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Abstract 

In this paper we give an exact analytical expression for the number of span- 
ning trees of an infinite family of outerplanar, small- world and self-similar 
graphs. This number is an important graph invariant related to different 
topological and dynamic properties of the graph, such as its reliability, syn- 
chronization capability and diffusion properties. The calculation of the num- 
ber of spanning trees is a demanding and difficult task, in particular for large 
graphs, and thus there is much interest in obtaining closed expressions for 
relevant infinite graph families. We have also calculated the spanning tree 
entropy of the graphs which we have compared with those for graphs with 
the same average degree. 
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1. Introduction 



Finding the number of spanning trees of graph is an old problem relevant 
in areas as diverse as mathematics, chemistry, physics, and computer science. 
This graph invariant is a parameter related, for example, to the reliability 
of a network 0| its synchronization 14] and the study of random walks 



[jj. The number of spanning trees of any finite graph can be computed 
from the well known Kirchhoff's matrix-tree theorem as the product of all 
nonzero eigenvalues of the Laplacian matrix of the graph [H] , however this is a 
demanding and difficult task, in particular for large graphs. Thus, there has 
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been much interest in finding other methods to produce exact expressions 
for the number of spanning trees of relevant graph families such as grids 
11 1, lattices 13] and Sierpinski gaskets Q, 

In this paper we give an exact analytical expression for the number of 
spanning trees of an infinite family of outerplanar, small-world and self- 
similar graphs which was introduced by two of the authors in This 
family also has clustering zero and an exponential degree distribution. The 
combination of modularity, small clustering coefficient, small- world prop- 
erties and exponential degree distribution can be found in some networks 
associated with several social and technological systems @, [Hj. Moreover, 
outer planarity is a relevant property as it is known that many algorithms 
that are NP-complete for general graphs perform polynomial in outerplanar 
graphs (l3 |. 



2. The M(t) graph family: construction and properties 

This family of graphs M(t) was introduced in 0] and can be constructed 
following a recursive- modular method, see Fig. [TJ 
For t = 0, M(0) has two vertices and one edge. 

For t = 1, M(l) is obtained from two graphs M(0) connected by two 
new edges. 

For t > 2, M(t) is obtained from two graphs M(t— 1) by connecting them 
with two new edges. In each M(i— 1) the two vertices chosen are adjacent 
with maximum degree and have also been used at step t— 1 to connect two 
M(i-2). 




M(0) M(\) M(2) M(3) M{f) 

Figure 1: Graphs M(t) produced at iterations t = 0, 1, 2, 3 and t. 

The number of vertices and edges of M{t) = (Vt,E t ) are, respectively, 
\V t \ = 2 t+1 and \E t \ = 3 • 2* - 2. 
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The graph M(t) is outerplanar, that is, the graph is planar and has an 
embedding in the plane such that all vertices lie always in the exterior face 
of the graph while the edges never cross. This family of graphs has the 
degrees of the endvertices positively correlated -the graphs are assortative- 
while their degree distribution follows an exponential law P CU m{k) = 2 2 ~ k . 
See 0, 0j for more details and properties. 

For large t, the diameter and average distance of M(t) scale logarithmi- 
cally with the order of the graph, meaning that the graph is small-world, 
although its clustering coefficient is zero. 

3. The number of spanning trees of M(t) 

In this section we make use of the construction method of M(t) described 
in the previous section to obtain a set of recursive equations for the number 
of spanning trees and spanning subgraphs of the graph, which then can be 
solved by induction to find the exact number of spanning trees of M(t) at 
any iteration step t. For this calculation, we adapt the decimation technique 
given in 0], which has been also considered to find the number of spanning 
trees of the Sierpinski gasket, the pseudofractal web, and several fractal 
lattices, see 0, [3, BJ] . 

As shown in the previous section, M{t) can be obtained from two M(t— 1) 
graphs interconnected by adding two new edges. As the four endvertices 
of these new edges are particularly relevant when counting the number of 
spanning trees of M(t), we denote them as at, bt, ct and dt and we will refer 
to them as hubs, see Fig. [TJ 

We recall here that a spanning subgraph of M(t) is a subgraph with the 
same vertex set as M(t) and a number of edges \E' t \ such that \E' t \ < \Et\. 
A spanning tree of M{t) is a spanning subgraph which is a tree and thus 
\E' t \ = \V t \-l. 

Let s(t) denote the total number of spanning trees of M{t) and let g(t) 
be the number of spanning subgraphs of M(t) which consist of two trees 
such that one of the hub vertices used to form M{t+1) belongs to one tree 
and the other is in the second tree (note that one of the trees might be just 
one hub). 

To find s(t), we classify the different cases that contribute to s(t) taking 
into account the connections between hub vertices of the two copies of M(t— 
1) that form M(t). 

We denote the existence of an edge between hubs a t , c t by {at, c t ) and we 
represent that there is no edge between them by (oj,q). We use the same 
notation for hubs bt and dt- 
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1. If (at,ct) and (bt,d t ) then there is no contribution to s(t). 

2. If (at, ct) and (bt, dt), the number of spanning trees that contribute to 
s(t) is 

s(t-l) 2 . (1) 

3. If (at, q) and (bt, dt), like in the previous case, the number of spanning 
trees is 

s(t-l) 2 . (2) 

4. If (at,ct) and (bt,dt). In this case the number of spanning trees con- 
tributing to s(t) is s(t— l)g(t— 1), i.e. the product of the number of 
spanning trees of the copy M(t — 1) containing hubs at and bt with 
g(t—l), which counts the number of spanning subgraphs of M(t—1) 
such that each contains the hubs q and dt but does not have a path 
connecting them because they belong to different trees of the spanning 
graph. We have also to consider the symmetrical case with respect to 
the sets {at, bt} and {ct, d t }. Thus, the total number of spanning trees 
that corresponds to this case is 

2s(t-l)g(t-l). (3) 

By adding ([T]), (J2J) and (|3J) we find that the total number of spanning 
trees of M(t) is 

s(t) = 2s(t-lf + 2s(t-l)g(t-l) (4) 

Now we need to obtain a recursive relation for g(t). A similar reasoning 
than above leads to the following equation 

g(t) = s(t-l) 2 + 2s(t-l)g(t-l). (5) 

Subtracting equation ([SJ from equation Q we obtain 

s(t) - g(t) = s(t-lf 

and thus 

g(t-l) = s(t-l) - s(t-2) 2 . (6) 
Introducing (jSJ) into ^ we have: 

s(t) = is(t-l) 2 - 2s(t-l)s(t-2) 2 with s(0) = 1 and s(l) = 4. 
To solve this equation, we rewrite it as 

s(t) 

s(P2) 



s(t-i) 2 im 



(7) 
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and, if we introduce 



q(t) = s(t)/s(t-l) 2 , t>2, 
?(1) = 4, 



equation (|7|) becomes 



q(t) = 4 



which can be solved to obtain 



g(t-l) 



q(t) = 2-V2 + 



, with t > 2, 



2^2 



1- (3 — 2^2)* " 
From equation (JSJ) we have: 

s(t) = q(t)s(t-l) 2 = q(t)q(t-l) 2 s(t-2) 4 = ■■■ 

= q(t)q(t-l) 2 q(t-2fq(t-3f ■ ■ ■ q(t-zf • • • ^lf" 1 

and thus, using (fTUl) . we can write: 

Theorem 3.1. The number of spanning trees of M(t), t > 1, is 

2V2 



'«) = n 



i=i 



2-V2+ 



1 - (3 - 2y/2)* 



(8) 



(9) 



(10) 



4. Spanning tree entropy of M(t) 

After having an exact expression for the number of spanning trees of 
M(t), we can calculate its spanning tree entropy, which is defined as in 0, 

lns(t) 



h = lim 



and we obtain 



Vt-»oo Vt 



h ~ 0.657. 



We can now compare this asymptotic value of the entropy of the span- 
ning trees of M(t) with those of other graph families with the same average 
degree: As an example, the value for the honeycomb lattice is 0.807 17J and 
the 4-8-8 (bathroom tile) and 3-12-12 lattices have entropy values 0.787 and 
0.721, respectively |l3l | while the spanning tree entropy for Hanoi graphs 
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is 0.677 [20(] . Thus, the asymptotic value of the entropy of the spanning 
trees of M(t) is the lowest known for graphs with average degree 3. This 
means that the number of spanning trees in M(t), although growing expo- 
nentially, do it at a lower rate than graphs with the same average degree. 
This result suggests that the family of graphs M (i) would be less reliable 
to a random removal of edges than the graphs mentioned above (as they 
have more spanning trees). However, the degree distribution also affects the 
reliability of a graph and it is well known that graphs with a scale-free dis- 
tribution of degrees are much more resilient than homogeneous graphs (like 
regular and exponential graphs, where all vertices are statistically identi- 
cal), see & Thus, the large inhomogeneity in the de gree distribution of 
the graphs M(t) might increase their robustness and it would be of interest 
to study, for these and other relevant graph families, the connections among 
spanning tree entropy, degree distribution and other topological parameters 
like degree correlation, in relation to their dynamic properties (reliability, 
synchronizability, etc.) 

5. Conclusion 

In this paper we find the number of spanning trees of the M(t) graph 
family by using a method, based on their self-similar structure, which allows 
us to obtain an analytical exact expression valid for any number of iterations. 
From the number of spanning trees we find their spanning tree entropy which 
is lower than in other graph families with the same average degree. 
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